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Abstract 



We re-examine neutrino oscillations in the early universe. Contrary to previous 
studies, we show that large neutrino asymmetries can arise due to oscillations 
between ordinary neutrinos and sterile neutrinos. This means that the Big Bang 
Nucleosynthesis (BBN) bounds on the mass and mixing of ordinary neutrinos with 
sterile neutrinos can be evaded. Also, it is possible that the neutrino asymmetries 
can be large (i.e. ~ 10%), and hence have a significant effect on BBN through 
nuclear reaction rates. 



There are several experimental indications that neutrinos have nonzero mass 
and oscillate [I], [| [|. It is possible that sterile neutrinos exist which oscillate with 
the known neutrinos. There are essentially two types of sterile neutrinos that can 
be envisaged. First, there are sterile states which either have no gauge interac- 
tions, or interactions which are much weaker than the usual weak interactions 
0]. Alternatively, it is possible to envisage neutrinos which do not have signifi- 
cant interactions with ordinary matter but do have significant interactions with 
themselves. An interesting example of the latter is given by mirror neutrinos|Q. 

However, for both sterile and mirror neutrinos there are apparently quite 
stringent bounds if they are required to be consistent with standard big bang 
cosmology. Assuming that the effective number of neutrino species present during 
Big Bang Nucleosynthesis (BBN) is bounded to be less than 3.4, then the mixing 
angle (8 ) and the squared mass difference (6m 2 ) for a sterile neutrino mixing 
with one of the known neutrinos is approximately bounded by 

iWKsin^o) 1 ' 8 ~ 10" 8 eV 2 , v = u e , 

|5m 2 |(sin 2 2fl ) L6 ~ KT 7 eV 2 , v = u^ T , (1) 

assuming that 5m 2 < and sin 2 26*o ~ 1CT 3 (for more precise bounds see e.g. 
figure 4 of Ref. f7j). These bounds arise by demanding that oscillations not bring 
the sterile neutrino into equilibrium with the known neutrinos. These bounds 
(along with even more stringent bounds for the case 5m 2 > or sin 2 2# ~ 1) 
would appear to exclude the region of parameter space required to explain the 
atmospheric neutrino anomaly in terms of — v s oscillation (|5m 2 | ~ 10~ 2 eV 2 , 
sin 2 26*o — 1), and would restrict the parameter space required to explain the 
solar neutrino deficit in terms of v e — v s oscillation. An important assumption in 
deriving the bounds of Eq.(l) was that the relic neutrino asymmetries could be 
neglected. The purpose of this paper is to re-examine this issue. We will show 
that for a significant region of parameter space large neutrino asymmetries can 
be generated due to neutrino oscillations. This has two important consequences. 
First, the bounds in Eq.(l) will be modified. Second, the asymmetries can be 
large enough to affect the neutron/proton ratio and hence also modify BBN. 
This is one plausible way of reconciling the possible disagreement of the BBN 
predictions with observations [§]. 

To define our notation, we first examine the ordinary neutrino (z/ Q , a = e, fi, r) 
oscillating with a sterile neutrino (v s ) in vacuum. Oscillations can occur if the 
weak eigenstate neutrino and sterile neutrino are each linear combinations 

v a = cos^qz/i + sin 6*0^2, 
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v s = — sin 6 i>i + cos 9 u 2 , (2) 

of mass eigenstates V\^. An ordinary neutrino of momentum p will then oscillate 
in vacuum after a time t with probability 

IM*)h)| 2 = sin 2 2fl sin 2 (-M , (3) 



where M 

2p _ 1 



(4) 



Sm 2 Aq 

However, in the early universe oscillations occur in a plasma. For v a — v s oscilla- 
tions in a plasma of temperature T, the matter and vacuum oscillation parameters 
are related by |L0 

• 2 9 n sin 2 26 

sm zv m - 1 _ 2zcQs29o + z 2i 

A 2 m = A 2 (l -2zcos26 + z 2 ), (5) 

where z = 2(p)(V a — V s )/5m 2 and (V a;S ) are the effective potentials due to the 
interactions of the neutrinos with matter ((p) ~ 3.15T). For a truly sterile 
neutrino Vg = 0. For a weak eigenstate neutrino u a , V a is given by O, O, III] 



^ Q = V2G F N, - ^!p) , (6) 



where Gf is the Fermi coupling constant, My/ is the W boson mass, A a is a 
numerical factor given by A e ~ 55 and A^ )T ~ 15.3 O, O. The quantity is 
given by 

L (a) = L Q + L Ve + L, M + L Ut + (±1 + 2x^)L e + - 2s w )L p - (7) 

where the plus sign refers to a = e and the minus sign refers to a = //, r. Also, 
a; m = sin 2 ^„, ~ 0.23 and L Q = (n Q -n s )/n 7 . 

We now study the possible creation of lepton number due to oscillations. This 
issue has been studied previously [12|, |13j where it was concluded that significant 



neutrino asymmetries could not be generated due to neutrino oscillations. That 
conclusion was based on studies covering a limited region of parameter space 
(in particular for very small values of \5m 2 \ ~ 10~ 5 eV 2 ). We will show below 
that for a large range of parameters it is possible to generate large neutrino 
asymmetries. Note that in Ref. |13| it was observed that the neutrino asymmetry 



could experience a brief period of exponential growth for appropriate parameters. 
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However, the particular parameters they chose led to a final asymmetry of only 
about 10~ 7 . The analysis below will demonstrate that much larger asymmetries 
(up to 10 _1 and even larger) can be generated for other interesting parameter 
choices. 

We will work analytically, and then check our work numerically using the 
density matrix formalism. Our analytic treatment assumes that the change in 
lepton number is dominated by collisions, which is true for temperatures greater 
than a few MeVfjHJ. We will also assume that the evolution of the neutrino 
ensemble follows the evolution of the state with average momentum. 

The change in lepton number due to collisions can be expressed in terms of 
reaction rates as follows [[y|: 

* (n "°- rag ° ) ~ -n Va T(u a — v s ) + n u Y(u s -> u a ) , g , 
+ n D T(v a -> v s ) - n D T(v s -> u a ) 

If the number densities of ordinary neutrinos are equal to the number densities 
of sterile neutrinos, then SL Ua /St — ► 0. This is because the rate Y[y a — > u s ) 
equals V(y B — > u a ) (and similarly for the antiparticle rates). However, if the 
number of sterile neutrinos is negligible to start with, and if they are not brought 
into equilibrium with the ordinary neutrinos, then SL Va /St is nonzero in general 
(provided Sm 2 < 0,008 26*0 > 0, as we shall demonstrate). We thus consider 
parameters satisfying the bounds Eq.(l), so that the sterile neutrinos are not 
brought into equilibrium with the ordinary neutrinos. 

Thus, assuming n Us , n Ps <C n Va , n 9a , the change in lepton number due to 
collisions can be expressed as, 



^ ^ \ [-I> Q - v s ) + T(u a - p.)] - [V(u a - u s ) + T(u a - p.)] , (9) 

where we have used n Va = 3n 7 /8 [recall L Ua = {n Va — n Va )/n^}. It will turn out 
that the second term (proportional to L Ua ) can be neglected for T ~ lMeV. 

Thus, in order to evaluate 8L Va / St, we need to evaluate the reaction rates. The 
rate T(v a — > v s ) is given by the interaction rate of ordinary neutrinos multiplied 
by the probability that the neutrino collapses to the sterile eigenfunction, i.e. 

Y{v a -> U s ) = {Pu a ^v s )colFv a (10) 

where T Ua = y a G 2 F T 5 (with y e ~ 4.0, y^ T ~ 2.9) are the reaction rates |J. The 
quantity (P Va -> Vs )coii is the probability that the neutrino u a collapses to the sterile 
state v s . It is given by 

(Pu a ^ s )coii = sin 2 2£ m (sin 2 -4%), (H) 

-Lose 
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where x is the distance between collisions. Note that (x) = Li nt = i/T Ua 
where Li nt is the mean distance between interactions. In the region where 
L mt » LH(IW), (sin 2 x/L(-)) - 1/2 and (sin 2 x/Z(™)> - 1/2 (ig and 
^osc being the oscillation lengths for neutrinos and antineutrinos, respectively). 
Thus, evaluating 5L Ua /St in the region where L mt > LH(LH), we find that@ 



^ ~ Ar Va (- sin 2 26 ma + sin 2 2^) . (12) 

Recall that sin 2 26 mct>m _ are defined in Eq.(5), where z is given by 

z = —a + b, for v = u a and 

z = a + b, for v — u a , (13) 

where 

a = -V^GpnjL^/Ao, b = -V2G F n 7 A a T 2 /(A M^). (14) 
Equation (12) then becomes 

8L Va 3T Ua sin 2 2# a(cos 26 - b) 



5t " 4[1 - 2 cos 29 (-a + b) + (a - b) 2 } [1-2 cos 26 (a + b) + (a + b) 2 ] 



(15) 



This contribution to 5L Ua /5t is larger than the term which we neglected [the 
second term of Eq.(9)] provided that a/(l + a 2 ) ^> AL Ua /3 cos 26*0, i.e. T ~ 
O.4(|5m 2 |/cos20 o eV 2 ) 1/4 - Thus, for T ~ lMeV the second term of Eq.(9) can 
be approximately neglected and Eq.(15) should be correct. Note that for tem- 
peratures of order 1 MeV and below, the effect of the second term in Eq.(9) is to 
reduce \L„ a \. 

We will assume that 5m 2 < (and hence a,b > 0) and cos 2^o > 00. 



Observe that for b > cos2# , = is a stable fixed point. To see this, note 
that if > then SL^/St < while for < 0, SL^/St > 0. The critical 
observation is that when b < cos2#o, = becomes an unstable fixed point 
(i.e. if i>) > then SL^/St > while for < 0, SL^/St < 0). Since b is 
proportional to T 6 , at some point during the evolution of the universe b becomes 
less than cos 29 and L (a) = becomes unstable. Observe that for a <C cos26*o, 
8L Va /8t is proportional to L Ua leading to a rapid exponential growth of L Ua . 
Calculating the temperature when b = cos2#o, we find 

, /cos 2d 1 5m 2 1 \ 1/6 . . 

T c ~ 12.9(16.2) f -£L j MeV, (16) 
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for electron (muon/tau) neutrinos. 

In order to calculate the amount of L Va produced, it turns out that it is 
more convenient to use the variable a. This is because Eq.(15) has simple forms 
depending on whether a < 1 or a > 1. The quantity 5a/ St, unlike 5L Va /5t, has 
a contribution from the expansion of the universe. Calculating 5a/ St, we find, 

5a da, da. dT da da dL Vty 

"777 osc "777 7777^ 



St dt lexp dt losc dt dT ' dL Ua dt 
-5.5T 3 4a 2V2G F n~ / 5L h 



(17) 



M P T A St ' 

where we have used the time-temperature relation t ~ Mp/llT 2 which is relevent 
for low temperatures (i.e. 2m e < T < m M ). We have also used the result that a 
is proportional to T 4 (so that da/dT = Aa/T) and that a is linear in L Ua . Thus, 
using Eq.(15) we find that 



5a 21T 2 a 



St ~ M f 



XT 7 (cos29 ~b) 



[(cos 26 + a - b) 2 + sin 2 26 ] [(cos 26 - a - b) 2 + sin 2 26 ] 

(18) 

where A is a constant given by 

-0.1ly a G F sm 2 26 M P 
~ 5^ 2 • (19) 

Note that A > since we are assuming that 5m 2 < 0. Observe also that 5a /St > 
for < a < a c (where a c is a number greater than cos26* which we will define 
precisely below), provided that T 7 cos26> > 1/A for b cos 29 (for b ~ cos2# , 
i.e. at resonance, the above behaviour holds for even smaller values of T). For 
a ^> cos 20q, 5a I St < 0. This means that the parameter a evolves towards a 
non-zero value: a — > a c ~ cos2# . To calculate a c we solve the equation 

£ = 0, (20) 

to find that 



, K + VK 2 - AC 



where 



< = —^~ 2 (21) 

K = 2(cos 26 - 6) 2 - 2 sin 2 26 , 
C= (1 -2cos26 b + b 2 ) 2 - XT 7 (cos 29 - b). (22) 
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Recall that L Va is related to a through Eq.(14). Hence, the result that a evolves 
to a c is equivalent to the statement that L Ua evolves to a non-zero value L c v ^ 
where L c v is given by 

K = t-° ac ■ (23) 

To summarise the situation: L Ua evolves to a very small value (so that — > 0) 
until the temperature is such that b = cos2# - At this point = becomes an 
unstable fixed point and L Ua exponentially increases until some point a c ~ cos 26 q. 
Strictly this behaviour will not occur for all values of 8m , sin 26 . For sin 2 26 
small enough, the evolution of a will be such that it will never reach the value a c . 
Nevertheless, there is a large range of parameters for which a evolves to a c . 
For later times, a c remains greater than cos 26q until 

M^E > max (^^j . (24) 

The maximum value of 8a\ osc /8t occurs at the resonance i.e. a ~ cos2# - Solving 
the above equation we find that 111 



For cos 26q ~ 1 then 



M P y a G 3 F cos 26 
\\8m 2 \/eV 2 )\ 117 



MeV. (26) 



Thus a c > cos 2^0 ~ 1 for T > T x . For later times (T < T x ), 8a\ exp /8t > 
max(8a\ osc /8t) and L Ua is essentially frozen. The relation a c ~ cos 28 and T = T x 
implies that the corresponding value for L Ua will be 

L Ua ^2xlO~ 2 (\8m 2 \/eV 2 )(T x )- i . (27) 

Recall that our analysis is really only valid for temperatures greater than about 
1 MeV. For \8m 2 \ < 30 eV 2 we should take T x to be of order lMeV. 

From Eq.(27) it is clear that we expect that large neutrino asymmetries can 
be generated. In fact, for \8m 2 \ ~ 5 eV 2 , L Ua ~ 0.1. We have checked our results 
numerically using the density matrix formalism (see e.g.fjHJ for a description of 
this formalism and for original references). Numerically integrating the equations 
describing the evolution of the density matrix we have found good agreement with 
the analytic results presented in this paper. In figure 1 we plot the evolution of L Va 
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taking, by way of example, Sm 2 = — 1 eV 2 , sin 2 29o — 10~ 8 . Also plotted in figure 
1 is the asymptotic value L c v [Eq.(23)]. Figure 1 shows the behaviour expected 
from our analytic analysis. The lepton number L Ve initially grows exponentially 
at the resonance b = cos 28 [which corresponds to 12.9 MeV for Sm 2 = —1 eV 2 , 
sin 2 26*o = 10~ 8 according to Eq.(16)]. After the exponential phase gets cut 
off, the neutrino asymmetry continues to grow more slowly as it approaches the 
asymptotic curve determined from a c . 

One important result is that the lepton number generated can be large enough 



to significantly modify the standard BBN scenario. Indeed if L Ue is positive |j20 



then it can reduce the neutron to proton ratio. (Such a reduction may be one 
way to get agreement with the data [§])■ 

Another important result of this work is that the bounds on ordinary neutrinos 
mixing with sterile neutrinos may be aleviated. We will study this issue in more 
detail and expand our analysis in a forthcoming article using the density matrix 
formalism. Below we sketch how this can happen. By way of example, assume 
that there is a sterile neutrino which mixes with the muon neutrino, with pa- 
rameters suggested by the atmospheric neutrino anomaly (i.e. \5m 2 \ ~ lO~ 2 eV 2 , 
sin 2 2^0 — 1@). Naively, this possibility is in conflict with BBN, because these 
parameters severely violate the BBN bounds of Eq.(l). However a small mixing 
between the tau neutrinos and the sterile neutrino may generate significant lep- 
ton number to strongly suppress the oscillations of the muon and sterile neutrino 



21|j . This would then make this scenario consistent with BBN. 

Consider a system comprised of z/ M , v T and v s . Note that for oscillations of 
tau to sterile neutrinos, the mass difference can be much bigger than the 10~ 2 eV 2 
for muon to sterile neutrinos. [<5m 2 s = m 2 s — m 2 ^ ~ — m 2 J. Denote 6 T (6 M ) 



as the value of the b parameter [defined in Eq.(14)] with 5m = 5m 2 s (5m 
Because of the larger mass difference and the larger value of cos20 TS , the point 
b T = cos 28 TS will be reached earlier for v T — v s oscillations than the corresponding 
point b^ = cos 26 ^ s for — v s oscillations. If the vacuum parameters for v T — v s 
mixing, 8 TS ,5m 2 s satisfy the BBN bounds of Eq.(l), then u s cannot be brought 
into equilibrium by v T — v s oscillation. Consequently our earlier analysis applies 
and we would expect significant generation of tau lepton number at this point. 
Note that the creation of L Vt could conceivably be compensated by — > v s 
oscillations. We assume that there is significant intergenerational mixing between 
Up and u T which we assume rapidly converts L Vr into so that L Ufi ~ L Ut . This 
assumption may not really be necessary (we will study this in more detail in our 
forthcoming article |p2"|). The condition that the creation of lepton number due 
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to v T — v s oscillations will not be erased by — u s oscillations is 

1^1 > 1^1, (28) 

for temperatures satisfing b T < cos 26 TS . Studying the above inequality away from 
resonance by using Eq.(15) (with b^/b T = 5m 2 s /5m 2 ,,), we find that it is sufficient 
to demand that 

sin 2 2^ s >sin 2 2^J^-f ) . (29) 



TS , 



If sin 2 29(j, s « 1 and Sm 2 s w 10 2 eV 2 then 

eV A 

sm 2 26 TS >lVr i —. (30) 

Thus, we arrive at a lower bound for sin 28 TS , which for m Vr < 30eV (as suggested 
by cosmology) is given by sin 

We also need to check that at the point where the lepton number begins 
to be generated [i.e. when b T = cos 29 TS , see Eq.(16)], the parameters are such 
that the v s have not already been brought into equilibrium. Demanding that 
\Y Va sin 2 26 m < H ~ 5.5T 2 /M P for T > T c [defined in Eq.(16)], we find: 



8m us I < cos 29 T 



1 5m 2 TS | 41 [cos 29 TS \ Sm 2 s |] i sin 2^ s ' 



(31) 



So, for |5m 2 s | ~ 10 2 eV 2 ,cos26 l rs ~ 1 and sin 26^ ~ 1, we get that m Vr ~ 
0.5 eV. 

Finally we need to check that the lepton number generated is sufficient to 
suppress the v^ — Vs oscillations. Demanding that interactions not bring the sterile 
neutrino into equilibrium with the muon neutrino, i.e. Y Vs ~ H ~ 5.5T 2 /Mp, we 
find in the case of large |21] that 

n A 

where T > Trf ec ~ 4AMeV (since we only need to require that the sterile neu- 
trinos not come into equilibrium before kinetic decoupling of the muon neutrinos 
occcur). Since from our earlier analysis, we expect a c ~ cos20 TS , we know that 
the generated by the u T — v s oscillations satisfies 

(r) > \5mU 

2 v / 2G F n 7 6.3T 
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Thus assuming that intergenerational mixing rapidly distributes the lepton num- 
ber so that ~ L^> ~ L^ T > then Eq.(32) is satisfied provided that 

\5m 2 TS \ Z 10|5m 2 J, (34) 

where the most stringent bound occurs when at the decoupling temperature i.e. 
T = Tjec — 4:AMeV. In other words, for the parameters satisfying the above 
equation there is sufficient lepton number generated to prevent the sterile neu- 
trino from coming into equilibrium (above T^ ec ) with the muon neutrino [despite 
having oscillation parameters in gross violation of the bounds Eq.(l)]. Assuming 
Sm 2 s ~ —m1 T and |5m 2 s | ~ 1CT 2 eV 2 (as suggested by the atmospheric neutrino 

anomaly) then Eq.(34) suggests that vtl Vt ~ 0.3 eV. This requirement is slightly 
less stringent then the 0.5 eV bound obtained earlier from Eq.(31). Thus we con- 
clude that muon neutrino oscillating into a sterile neutrino with maximal mixing 
and 5m 2 ~ 10~ 2 eV 2 as suggested by the atmospheric neutrino anomaly seems to 
be consistent with BBN provided that the tau neutrino oscillates into the sterile 
neutrinos with parameters bounded by 10~ 4 /m^ T < sin 2 28 TS ~ 3 x lO^/m*- 2 
[where we have used Eqs.(l) and (30)]. Note that this condition implies that 
m Vr > 1.5 eV. 

This result should also apply for mirror neutrinos which do not have signifi- 
cant interactions with ordinary particles but do have significant interactions with 
themselves 0. Mirror neutrinos are required to exist if there is an unbroken par- 
ity symmetry. However, naively this scenario appears to be in conflict with BBN 
because it predicts that the muon neutrino will be maximally mixed with a mir- 
ror neutrino. Thus it will violate the bounds Eq.(l) if it is to be relevant for the 
atmospheric neutrino anomaly. However, we have shown in this paper that this 
scenario is not in conflict with BBN provided that the sterile neutrino also mixes 
slightly with the tau neutrino. 

Thus we conclude that it is possible that large neutrino asymmetries can be 
generated by oscillations between ordinary and sterile neutrinos. This asymmetry 
can be very large (i.e. ~ 0.10) and may thus modify nucleosynthesis. This should 
provide one way to reconcile the present discrepancy between the BBN predictions 
and observations ||. Another consequence of a large neutrino asymmetry is that 
it may be possible to aleviate the BBN bounds on the mixing of ordinary and 
sterile neutrinos. 
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Figure Caption 

Electron lepton number (L Ue ) versus temperature (in MeV units) for v e — u s 
oscillations with parameters Sm 2 = — 1 eV 2 and sin 2 29 = 10~ 8 . The broken line 
is the theoretical prediction based on the numerical integration of the density 
matrix equations. The solid line is the asymptotic value L c Ue [Eq.(23)]. 
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